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Semileptouic B decays to DXiV {£ = e or /i) are selected by reconstructing D"£ and D^£ combi- 
nations from a sample of 230 million T(45) BB decays recorded with the BABAR detector at the 
PEP-II e+e^ collider at SLAC. A global fit to these samples in a 3-dimensional space of kinematic 
variables is used to determine the branching fractions B{B~ — » D^tP) = (2.34 ± 0.03 ± 0.13)% and 
B{B- -* D*°ep) = (5.40 ± 0.02 ± 0.21)% where the errors are statistical and systematic, respec- 
tively. The fit also determines form factor parameters in a HQET-based parameterization, resulting 
in pi, = 1.20 ± 0.04 ± 0.07 for B D(17 and = 1.22 ± 0.02 ± 0.07 for B D*eM. These values 
are used to obtain the product of the CKM matrix element |Vc6| times the form factor at the zero 
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recoil point for both B Dgp decays, G{l)\Vcb\ = (43.1 ± 0.8 ± 2.3) x 10"^, and for B -» D*£l7 

decays, J^{l)\Vcb\ = (35.9 ± 0.2 ± 1.2) x 10"^ 



PACS numbers; 13.25.Hw, 12.15.Hh, ll.30.Er 
I. INTRODUCTION 

The study of scmilcptonic decays of heavy quarks pro- 
vides the cleanest avenue for the determination of sev- 
eral elements of the Cabibbo-Kobayashi-Maskawa ma- 
trix [1], which arc fundamental parameters in the stan- 
dard model of particle physics. The coupling strength of 
the weak b c transition is proportional to \Vcb\, which 
has been measured in both inclusive semileptonic B de- 
cays [2] and in the exclusive transitions B — > DiU [3-6] 
and B — > D*iP [3, 6-10] {£ = e or fi and charge con- 
jugate modes are implied). The inclusive and exclusive 
determinations of \Vcb\ rely on different theoretical calcu- 
lations. The former employs a parton-levcl calculation of 
the decay rate organized in a double expansion in as and 
in inverse powers of mi,, the &-quark mass. The latter re- 
lics on a parameterization of the decay form factors using 
Heavy Quark Symmetry and a non-perturbative calcula- 
tion of the form factor normalization at the zero recoil 
(maximum squared momentum transfer) point. The the- 
oretical uncertainties in these two approaches are inde- 
pendent. The inclusive and exclusive experimental mea- 
surements use different techniques and have negligible 
statistical overlap, and thus have largely uncorrelated un- 
certainties. This independence makes the comparison of 
\Vcb\ from inclusive and exclusive decays a powerful test 
of our understanding of semileptonic decays. The latest 
determinations [11] differ by more than two standard de- 
viations (cr), and the inclusive determination is currently 
more than twice as precise as the exclusive determination. 
Improvements in the measurements of exclusive decays 
will strengthen this test. This is particularly true for 
the B — ^ DiJy decay, where the experimental uncertain- 
ties dominate the determination of \Vcb\- For the decay 
B^ — > D*~^iU, the experimental situation needs clarifi- 
cation, as existing measurements arc in poor agreement 
with each other [11]. Finally, precise measurements of 
semileptonic B decays to charm are needed to further 
improve determinations of \Vub\, where B D^*^£V de- 
cays are the principal background. 

Semileptonic b ^ c transitions result in the produc- 
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tion of a charm system that cascades down to the ground 
state Z?" or D+ mesons. Most previous analyses have fo- 
cused on reconstructing separately the exclusive decays 
B D*(V [3, 7-10] and B DiU [3-5]. The B D*iU 
analyses involve reconstruction of the soft transition pion 
from the decay D* Dtt, which is at the limit of detec- 
tor acceptance; determination of the reconstruction ef- 
ficiency for these pions introduces significant systematic 
uncertainty. Studies of the exclusive decay B — > DiU 
suffer from large feed-down background from B D*iu 
decays where the transition pion is undetected. 

In this analysis we reconstruct D^t and D^i pairs and 
use a global fit to their kinematic properties to deter- 
mine the branching fractions and form factor parame- 
ters of the dominant semileptonic decays B — > Diu and 
B D*il7. The reconstructed D£ samples contain, by 
design, the feed-down from all the higher mass states 
(apart from decays of the type B — > DfXiu [12]). Kine- 
matic restrictions are imposed to reduce the contribution 
of backgrounds from semileptonic decays to final state 
hadronic systems more massive than D* and from other 
sources of D£ combinations. Distributions from selected 
events are binned in the 3-dimensional space described 
below. The electron and muon samples are input into 
separate fits, in which isospin symmetry is assumed for 
the semileptonic decay rates. Semileptonic decays are 
produced via a spectator diagram in which the heavy 
quark decays independently; strong interaction correc- 
tions to this process conserve isospin. As a result, we 
constrain semileptonic decay rates for B~ and B° to be 
equal, e.g., r{B- DH-p) = T{B° D+l'V). This 
substantially reduces statistical uncertainties on the fit- 
ted parameters. Systematic uncertainties associated with 
the modeling of the signal and background processes, 
the detector response, and uncertainties on input pa- 
rameters are determined, along with their correlations 
between the electron and muon samples. The fitted re- 
sults are then combined using the full covariancc matrix 
of statistical and systematic errors. For both B Div 
and B — > D*iT' decays, the fitted branching fractions and 
form factor parameters are used to determine the prod- 
ucts CJ(l)|ycfc| and JF(l))ycfc). These measurements, along 
with theoretical input on the form factor normalizations 
Q{\) and .^-"(1) at the zero recoil point, allow determina- 
tions of \Vcb\- 

The approach taken in this study has some similar- 
ity to that of Ref. [6], where the branching fractions 
for B — > Div and B ^ D*iu are measured simulta- 
neously. However, Ref. [6] reconstructs semileptonic B 
decays in events in which the second B meson is fully 
reconstructed. That approach allows the use of the miss- 
ing mass squared as a powerful discriminant. This analy- 
sis provides modest discrimination between the different 
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semileptonic decays on an cvcnt-by-cvcnt basis, but re- 
sults in a much larger statistical sample and enables the 
measurement of form factor parameters. 

The remaining sections of this paper arc organized as 
follows. In Sec. II we describe the BaBAR. detector and 
the samples of BaBAR data and simulated events used 
in the analysis. The event selection and the distribu- 
tions that are input to the global fit are discussed in 
Sec. III. Wc give the parameterization of the form fac- 
tors of S ^ decays and the modeling of semilep- 
tonic B decays to £>(*V and D^*^TnT states in Sec. IV. 
The global fit strategy and results are given in Sec. V, 
and the evaluation of systematic uncertainties is detailed 
in Sec. VI. Sec. VII presents the determination of \Vcb\ 
from the fitted results. The final section (VIII) discusses 
the results and provides averages with previous BaBAR 
measurements. 



II. THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with 
the BaBar detector at the PEP-II storage ring between 
1999 and 2004. PEP-II is an asymmetric collider; the 
center-of-mass of the colliding e^e^ moves with veloc- 
ity P = 0.49 along the beam axis in the laboratory rest 
frame. The data collected at energies near the peak of 
the T{AS) resonance (on-peak) correspond to 207 fb~^ 
or 230 million BB decays. Data collected just below BB 
threshold (off-peak), corresponding to 21.5 fb~^, are used 
to subtract the e+e" qq {q = u,d, s,c) background un- 
der the T{4S) resonance. 

The BaBAR detector is described in detail else- 
where [13]. It consists of a silicon vertex tracker (SVT), 
a drift chamber (DCH) , a detector of internally reflected 
Cherenkov light (DIRC), an electromagnetic calorimeter 
(EMC) and an instrumented flux return (IFR). The SVT 
and DCH operate in an axial magnetic field of 1.5 T and 
provide measurements of the positions and momenta of 
charged particles, as well as of their ionization energy 
loss {dE/dx). Energy and shower shape measurements 
for photons and electrons are provided by the EMC. The 
DIRC measures the angle of Cherenkov photons emitted 
by charged particles traversing the fused silica radiator 
bars. Charged particles that traverse the EMC and show- 
ering hadrons are measured in the IFR as they penetrate 
successive layers of the return yoke of the magnet. 

Simulated events used in the analysis are generated us- 
ing the EVTGEN [14] program, and the generated parti- 
cles are propagated through a model of the BABAR detec- 
tor with the GEANT4 [15] program and reconstructed 
using the same algorithms used on BaBar data. The 
form factor parameterization [16] used in the simulation 
for B — > D*ti' decays is based on Heavy Quark Effec- 
tive Theorj^(HQET) [17],_while the ISGW2 model [18] 
is used for B — > Div and B — > D**iP decays, where D** 
is one of the four P-wave charm mesons as described 
in Sec. IIIB. These are subsequently reweighted to the 



forms given in Sec. IV. For non-resonant B D^*'^T:tU 
decays, the Goity-Roberts model [19] is used. In order 
to saturate the inclusive semileptonic h ctP decay rate 
we include a contribution from B D^*^Tnr£T' decays; a 
variety of models are considered for this purpose. The 
branching fractions for B and charm decays in the simu- 
lation are rescaled to the values in Ref. [11]. In addition, 
the momentum spectra for and D+ from B DX 
and B DX decays are adjusted to agree with the 
corresponding measured spectra from Ref. [20] . This ad- 
justment is done only for background processes. 

The simulation of the detector response provided by 
the GEANT4-based program is further adjusted by com- 
paring with BaBar data control samples. In particular, 
the efficiency of charged track reconstruction is modi- 
fied by 1-2%, depending on momenta and event multi- 
plicity, based on studies of multi-hadron events and 1- 
versus-3 prong e+e~ t^t~ events. The efficiencies 
and misidentification probabilities of the particle identifi- 
cation (PID) algorithms used to select pions, kaons, elec- 
trons and muons (see Sec. HI) are adjusted based on stud- 
ies of samples of e+e~ — > e+e~7 and e+e~ — > /i+/Lt~7, 
and several samples reconstructed without particle iden- 
tification: l-versus-3 prong e+e~ t^t~ events, 
n+TT-, D*+ D°n+ {K-tt+)'jt+ and A pw'. 



III. EVENT SELECTION 

A. Preselection of D£ candidates 

We select multi-hadron events by requiring at least 
three good-quality charged tracks, a total reconstructed 
energy in the event exceeding 4.5 GeV, the second nor- 
malized Fox-Wolfram moment [21] i?2 < 0.5, and the 
distance between the interaction point and the primary 
vertex of the B decay to be less than 0.5 cm (6.0 cm) 
in the direction transverse (parallel) to the beam line. 
In these events an identified electron or muon candidate 
must be present, along with a candidate D meson de- 
cay. Candidate electrons are identified using a likelihood 
ratio based on the shower shape in the EMC, dE/dx in 
the tracking detectors, the Cherenkov angle and the ratio 
of EMC energy to track momentum. The electron iden- 
tification efficiency is 94% within the acceptance of the 
calorimeter, and the pion misidentification rate is 0.1%. 
Muon candidates are identified using a neural network 
that takes input information from the tracking detectors, 
EMC, and IFR. The muon identification efficiency rises 
with momentum to reach a plateau of 70% for laboratory 
momenta above 1.4 GeV/c, and the pion misidentification 
rate is 3%. 

Kaon candidates are required to satisfy particle iden- 
tification criteria based on the dE/dx measured in the 
tracking detectors and the Cherenkov angle measured 
in the DIRC. Each kaon candidate is combined with 

one or two charged tracks of opposite sign to form a 
D° K-TT+ or D+ K-'K+Tr+ candidate. Those 
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combinations with invariant masses in the range 1.840 < 
niKTr < 1-888 GeV/c^ are considered as 1?° candidates 
and those in the range 1.845 < mK-mr < 1.893 GeV/c^ 
as candidates, respectively. Combinations in the 
"sideband" mass regions 1.816 < mxTT < 1-840 GeV/c^ 
and 1.888 < ttik^ < 1.912 GeV/c^ (1.821 < m^^^ < 
1.845 GeV/c2 and 1.893 < m^TTTr < 1-917 GeV/c^) are 
used to estimate the combinatorial background. 

The charge of the kaon candidate is required to have 
the same sign as that of the candidate lepton. Each D- 
lepton combination in an event is fitted to both B — > 
D£ and D K~tt'^{tt^) vertices using the algorithm 
described in Ref. [22] . The fit probabilities are required to 
exceed 0.01 for the B D^l and B — > _D+£ vertices and 
0.001 for the and decay vertices. We require the 
absolute value of the cosine of the angle between the Dl 
momentum vector and the thrust axis of the remaining 
particles in the event to be smaller than 0.92 to further 
reduce background, most of which comes from e^e~ — *■ 
qq {q = u, d, s, c) events. 

The signal yields are determined by subtracting the es- 
timated combinatorial background from the number of D 
candidates in the peak region. The combinatorial back- 
ground is estimated using the number of candidates in 
the D mass sideband regions scaled by the ratio of the 
widths of the signal and sideband regions. This is equiva- 
lent to assuming a linear dependence of the combinatorial 
background on invariant mass. The change in the yields 
is negligible when using other assumptions for the back- 
ground shape. Candidates from e+e^ qq events are 
statistically removed from the data sample by subtract- 
ing the distribution of candidates observed in the data 
collected at energies below BB threshold (off-peak), after 
scaling these data by the factor rc = (-ConSoff) / ('CoffSon) 
to account for the difference in luminosity and the depen- 
dence of the annihilation cross-section on energy. The se- 
lection criteria listed above were determined using sim- 
ulated BB events and off-peak data to roughly maxi- 
mize the statistical significance of the Dl signal yields 
in e^e^ — > BB events. They have an overall efficiency 
of 80% (76%) for B D^X^v {B D+X^v) decays 
with p| , the lepton momentum magnitude in the center- 
of-mass (CM) frame, in the range 0.8-2.8 GeV/c. 

The invariant mass distributions for the Z?" and 
candidates, after off-peak subtraction, are shown in Fig. 1 
for two kinematic subsets representing regions with good 
and poor signal-to-background ratios. The small dif- 
ferences in peak position and combinatorial background 
level have negligible impact on the analysis due to the 
sideband subtraction described above and the wide sig- 
nal window. 

The D^l and candidates are binned in three kine- 
matic variables: 

• p*]j, the D momentum in the CM frame; 

• the lepton momentum in the CM frame; 



the cosine of the angle between the B and D£ mo- 
mentum vectors in the CM frame under the as- 
sumption that the B decayed to D£v. If the D£ 
pair is not from a B —> D£v decay, | cosOs-Del can 
exceed unity. The B energy and momentum are not 
measured event-by-event; they are calculated from 
the CM energy determined by the PEP-II beams 

as E*Q = and = \J E*^'^ - m\, where uib 

is the B^ meson mass. The energy, momentum and 
invariant mass corresponding to the sum of the D 
and lepton four vectors in the CM frame are de- 
noted E*jj^, p*jj^ and moi, respectively. 

The binning in these three variables is discussed in 
Sec. inc. 



6000^ 




m(K3i) (GeV/c^) m(K;c;i) (GeV/c^) 



FIG. 1: The invariant mass distributions (data points) for 
selected candidates. Scaled off-peak data have been sub- 
tracted to remove contributions from e'*"e~ — > qq annihila- 
tion. Plots (a,c) show K~tt'^ combinations and (b,d) show 
K'tv'^-k'^ combinations. In each case the D£ candidates are 
required to satisfy —2.0 < cosOb-di < l-l- The further 
kinematic requirements are 1.6 < pi < 1.8 GeV/c, 1.6 < 
ph < 2.0 GeV/c for plots (a,b) and 2.0 < Pe < 2.35 GeV/c, 
0.8 < pj) < 1.2 GeV/c for plots (c,d). The histograms show 
the contribution from simulated BB events scaled to the data 
luminosity. The arrows indicate the boundaries between sig- 
nal and sideband regions. 



B. Sources of D£ candidates 

There are several sources of D£ candidates that survive 
the ZJ-mass sideband and off-peak subtractions. In both 
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the and samples we group them as follows {B 
represents both B~ and B^): 

(i) B -> Dei7 

(ii) B ^ D*ev 

(iii) B — > D^*^{mr)£l7, which includes 

• The P-wave Z?** charm mesons. In the frame- 
work of HQET, the P-wave charm mesons are 
categorized by the angular momentum of the 
light constituent, jg, namely jf — l/2~ dou- 
blet and D[ and j[ = 3/2- doublet Di 
and [23]. 

• Non-resonant B D^*^tt^. 

• Decays of the type B D'^*^'K'KtP; the model- 
ing of these is discussed in Sec. IV D. 

(iv) Background from BB events in which the lepton 
and D candidates do not arise from a single semilep- 
tonic B decay. These include (in order of impor- 
tance) 

• Direct leptons from B X£v decays com- 
bined with a D from the decay of the other B 
meson in the event. Roughly one third of this 
background comes from events in which B^B^ 
mixing results in the decay of two B^ mesons. 
Most of the remaining contribution comes from 
CKM-suppressed B DX transitions. 

• Uncorrelated cascade decays. In this case the 
lepton mostly comes from the decay of an anti- 
charm meson produced in the B decay and the 
D arises from the decay of the other B meson 
in the event. 

• Correlated cascade decays, in which the lepton 
and D candidates come from the same parent 
B meson. These are mainly B — DD{X) and 
B D{X)tv decays, with the lepton coming 
from the decay of an anti-charm meson or tau. 

• Mis-identified lepton background. The prob- 
ability of a hadron being misidentified as a 
lepton is negligible for electrons but not for 
muons. 

As mentioned previously, the same decay widths are im- 
posed for the semileptonic transitions of B^ and B^ . For 
the background processes (source iv) no such requirement 
is imposed. 



C. Kinematic restrictions 

Despite the use of the best available information for 
calculating the background and B — > D'^*\n'K)£v dis- 
tributions, these components suffer from significant un- 
certainties. We therefore restrict the kinematic range 
of the variables used in the fit to reduce the impact 



of these uncertainties while preserving sensitivity to the 
B Dtp and B 0*117 branching fractions and form 
factor parameters. We require —2 < co&Ob-di < 1-1 
and place restrictions on p'^ and rejecting regions 
where the signal decays are not dominant. This re- 
sults in the ranges 1.2 GeV/c < p| < 2.35 GeV/c and 
0.8 GeV/c <ph< 2.25 GeV/c. The yield within this re- 
gion is 4.79 X 10^ (2.95 x 10^) candidates in the D"i [D+tj 
sample with a statistical uncertainty of 0.26% (0.66%). 



TABLE I: Definition of bins used for kinematic variables. 



Quantity ^ bins 



Bin edges 



cosOb-di 3 -2.0, -1.0, 0.0, 1.1 

p'l (GeV/c) 10 1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.1,2.35 

Pd (GeV/c) 8 0.8,1.1,1.35,1.5,1.65,1.8,1.95,2.1,2.25 

The data are binned finely enough to have good sen- 
sitivity to the fit parameters while maintaining adequate 
statistics per bin. Table I gives the binning used in the 
fit. We avoid setting a bin edge at cos 9 b-D£ = 1 to re- 
duce our sensitivity to the modeling of the resolution in 
this variable, since the B — > DiU decay distribution has 
a sharp cut-off at this point. 
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FIG. 2: (Color online) Distribution of pjj vs. for _D°e can- 
didates after sideband subtraction. The shaded boxes have 
area proportional to the number of entries. The plots show 
simulated candidates for (a) B — > DeV, (b) B — > D*eV and 
(c) other (sources iii and iv combined) , and for data after off- 
peak subtraction (d). The binning given in Table I is used 
and only candidates that satisfy 0.0 < cos8b-di < 1-1 are 
plotted. 



Two-dimensional projections of the signal, background 
and data distributions for the D'^e sample are shown in 
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FIG. 3: (Color online) Projections onto individual kinematic variables of the data after off-peak subtraction and the results of 
the fit : (a,d) lepton and (b,e) D momentum in the CM frame, and (c,f) cos9b-di- The points show data for accepted D^e 
(a,b,c) and D'^e (d,e,f) candidates, and the histograms show the individual fit components (from top to bottom): B — + DeV, 
B D*eV, B D^*\mT)eV and other BB background. The ratio of data to the sum of the fitted yields is shown below each 
plot. 



Fig. 2 to illustrate the separation power in these vari- 
ables. The distributions for the D'^fi sample (not shown) 
are similar. The one-dimensional projections of the De 
and Dfi samples are shown in Figs. 3 and 4. The differ- 
ence in the size of the B D*£v components in D'^l and 
distributions is due to the fact that D*^ does not 
decay io . 



IV. MODELING OF SEMILEPTONIC B 
DECAYS 

In our fully simulated event samples B DtU and 
B D**ii7 decays were generated using the ISGW2 
model [18]. For B D*ep decays, an HQET model 
was used with a linear form factor parameterization. We 
re- weight all these decays using the formulae given in the 



following subsections. The histograms in Figs. 3 and 4 
are re- weighted. 

A. B Dtp decays 

The differential decay rate is given by [17] 

dw 487r3 ^ ' 

x[{\ + r)h+{w)-(l-r)hPw)Y, (1) 

where Gp is the Fermi constant, h^{w) and h^(w) are 
the form factors, r = mo/mB is the mass ratio and ms 
and mo are the B and D meson masses, respectively. 
The velocity transfer w is defined as 

w = vb-vd, (2) 
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FIG. 4: (Color online) Projections onto individual kinematic variables of the data after off-peak subtraction and the results of 
the fit : (a,d) lepton and (b,e) D momentum in the CM frame, and (c,f) cos^s-of- The points show data for accepted D'^/i 
(a,b,c) and (d,e,f) candidates, and the histograms show the individual fit components (from top to bottom): B — > D/i77, 

B — > D* ^u, B D(*)(n7r) fiu and other BB background. The ratio of data to the sum of the fitted yields is shown below 
each plot. 



where vb and vd are the 4- velocities of the B and D 
mesons, respectively. In the B rest frame w corresponds 
to the Lorentz boost of the D meson. In the HQET 
model, the form factors are given by [16] 



[1 - 8plz + (51p2 



D 



10)^2 



84)z3 



(3) 



and 



h_{w)=0, (4) 

where z = {\/w + 1 — ^/2)/{^/w -\- 1 -I- v^) and pj^ and 
CJ(1) are, respectively, the form factor slope and normal- 
ization at w = 1. 

The above formulae neglect the lepton mass mi. Muon 
mass effects need to be included to achieve precision at 
the few percent level on the form factor parameters. Al- 
lowing for non-zero lepton mass introduces additional 



terms in the phase space and form factor expressions [24] 
that can be included by multiplying the decay rate for- 
mula by the following factor: 



Wd = 
where 
K.d{w) = 



1 + — 2rw rn?p. 



1 + ]Cd{w)- 



(5) 



1 -r 
1 + r 



w + 1 

w — 1 



1 



2(1 + r2 - 2rw) ' 



(6) 



B. B D'gu decays 

We need three additional kinematic variables to de- 
scribe this decay. A common choice is 6i, Oy and x. 
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shown in Fig. 5, and defined as 

• 9i : the angle between the lepton and the direction 
opposite the B meson in the W rest frame. 

• 0y : the angle between the D meson and the direc- 
tion opposite the B meson in the D* rest frame. 

• X ■ the azimuthal angle between the planes formed 
by the W-i and D* — D systems in the B rest frame. 




FIG. 5: Definition of the three angles 9i, 9v and x for the 
decay B D'H/. 



The form used for hAi{w) is [16] 
hA,{w) = ^(1) X 

[l - Spl.z + {53pl. - 15)z2 - (231p^. - 9l)z^] , (11) 

where and .F(l) are, respectively, the form factor 
slope and normalization ai w — 1. 

Non-zero lepton mass is accounted for by multiplying 
the decay rate formula by the factor 



where 
JCd' (w) 



1 + r*^ — 2r*w m?j 



2 r 



i + /Cz,.H^ 



(12) 



3 Hi 

2 Hi + Hi 



Hi 



1 



2(l + r*2 - 2r*w)' 
(13) 

Here, Ht is expressed, using another form factor ratio 
Rsiw), by 



The differential decay rate is given by [17] 

dr{B -> D*£i7) 
dw dcosOy dcosOf dx 
3Ge 



4(47r)' 



iKbl^msm^. y/w^ - 1(1 + r*2 - 2r*w) x 



[(l-cos6'f)^sin^6'v|7?+(w)|^ 
+ {1 + cosei^fsm^ev\H_{w)\^ 
+4sm^eecos'^ev\Ho{w)\^ 

— 4sin6'f(l — cos9i)sm6v cosdv cosxH-^-{w) Ho{w) 
+Asin9i{l + cos9e)sm9vcos9vcosxH^{w)Ho{w) 
-2sm^9eam^9vcos2xH+ {w)H^ (w)] , (7) 

where Hi{w) are form factors, r* = rriDt/mB and mo* 
is the D* meson mass. The Hi{w) are usually written 
in terms of one form factor Ha-^ (w) and two form factor 
ratios, Ri{w) and R2{w), as follows: 

R*(l - r*'^)(w + 1) 



2vr 



2r*w 



where R* = {2^mBmD') / {itib + m/j*) and 



H±{w) = ^ [l^J R,{nj) 



HoH - l + i^(l-i?2(z«)) 



w + 1 



(9) 



The form factor ratios have a modest dependence on w, 
estimated [16] as 

Ri{w) = i?i -0.12(w- l)-f-0.05(w- 1)2, 
R2iw) = i?2+0.11(u)-l)-0.06(w-l)^ (10) 



Htiw) 



1 _ r* 



r* — w 1 + r*2 — 2r*w 

1 H -^^13 — r~r^ — R2 



w + 1 



r*(w -I- 1) 



(14) 

We take Rsiw) — 1; this approximation has a negligible 
impact on our fit results. 



C. B ^ D'-*\£V decays 

The four P-wave D** states have been measured in 
semileptonic decays [25-27]. The decays B — > D**^^ are 
modeled following an HQET-inspired form factor param- 
eterization given in Ref. [23] . Detailed formulae are given 
in Appendix A. We use the approximation Bi of this 
model for our main fit and use the approximation B2 to 
evaluate the uncertainty due to the approximation. The 
slope of the form factors versus w is parameterized by f', 
which we set to —1.5 and vary between —1.0 and —2.0 
to study systematic uncertainties (Table II). 

To parameterize the B D^*^'Kti> decay branching 
fractions we define five branching fraction ratios: 



, _ B{B 

JDi/Di — 



JDtt/D-^ = 



ID', 



B'^^^iB- D*+7r-iV) 

B{B- Dfeu) 



(15) 



(16) 



(17) 



^ B{B- D*°eT7) + B^^jB- DTveu) 
Sdid^/d,di = 6(5- ^ D'itV) + B{B- ^ Df£u) ' 

(18) 
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JD[D'7t/DiD* = 



B{B~ Dleu) + B{B- D*°ev) ' 

(19) 

where NR stands for "non-resonant" decays, which are 
assumed to be isospin- invariant. The quantity fo'/Di is 
the ratio between two narrow states, Jd-k/d^ (/d-tt/d;) is 
between two broad states decaying to D-k {D*tt) and the 
other two ratios are between broad and narrow states. 
With these definitions the branching fractions for indi- 
vidual modes can be related to the total branching frac- 
tion B{B- £)(*)7r^) = B{B- Dntv) + B{B- 
D*tt£T'). We combine a new measurement [6] with the 
world average [11] to determine the value given in Ta- 
ble II. 

To estimate the branching fraction ratios, we average 
several measurements [25-28] to find 



B{B- 




= 0.0042 ± 0.0004, 




B{B- - 




= 0.0031 ± 0.0005, 


(20) 


B{B- - 




= 0.0022 ±0.0014, 




B{B- - 




= 0.0048 ± 0.0008. 





The sum of the D** branching fractions saturates 
B{B~ D'^*^'k£v) which implies that the non-resonant 
branching fractions are small. We use 



B^'^B- 
B^^{B- 



^ Dn(J7) = 0.0015 ±0.0015 
D*7r£l7) = 0.00045 ±0.00045. (21) 



From these numbers the branching fraction ratios are cal- 
culated and listed in Table II. These quantities are taken 
as independent when evaluating systematic uncertainties. 



D. B ->■ D'^*'>niT£P decays 

Recent measurements [6, 11] indicate that the inclu- 
sive B — > XciP branching fraction is not saturated by 
the sum of the B -> Dgp, B D*tU and B D'-*\eu 
branching fractions. In order to fill the gap, we include 
B — !■ D'^*^'K'Kti' decays in our fit. We assume the branch- 
ing fraction of these decays, given in Table II, is equal 
to this missing contribution to the inclusive branching 
fraction [11]. 

The B D(*)7r7r£F decays are modeled as a combina- 
tion of four resonances : pseudo-scalar {Xc) and vector 
{X*) states just above D*TnT threshold, and a heavier 
pair of pseudo-scalar (Fc) and vector {¥*) states just 
above D*p threshold, as listed in Table III. Each state 
is assumed to be produced with equal rate in semilep- 
tonic B decays and each is assumed to decay with equal 
branching fraction to Dtttt and 0*1^, conserving isospin. 
These assumptions are varied in assessing systematic un- 
certainties. 



V. GLOBAL FIT 

The binned distributions of D^i and D~^£ candidates 
in the variables p^, and cos 9 b-d£ arfJ fitted with 
the sum of distributions for the signal and background 
sources listed in Sec. IIIB. The expected shape of the 
individual components is based on simulation, and the 
fit adjusts the normalization of each component to min- 
imize the global chi-squared: 



bin i 



(22) 

where the index i sums over bins of the D^i and D'^l 
distributions and j sums over individual simulated com- 
ponents. The coefficients Cj depend on a, the set of free 
parameters determined by minimizing . For example, 
for the B~ ^ D'^tu component the coefficient Cj is given 
by the ratio of the fitted B{B~ — > D^tu) branching frac- 
tion to the value used in generating the corresponding 
distribution. The number of candidates in the data col- 
lected on (below) the T(4S') peak in bin i is denoted A^°° 
(7V°*) and Mij is the number of simulated events in bin 
i from source j. The may depend on a as explained 
below. The statistical uncertainties, after D mass side- 
band subtraction, are given by the ai for the data and 
the aij for the different Monte Carlo samples. The factor 
Tj is the ratio of the on-peak luminosity to the effective 
luminosity of the appropriate Monte Carlo sample. Only 
those bins in which the number of entries expected from 
the simulation exceeds 10 are used in the sum. 

For the B DtP and B D*£T' signal components 
we fit for both the branching fractions and for form-factor 
parameters. To facilitate this, we split these components 
into sub-components, one corresponding to each unique 
combination of the parameters a in the expression for 
the decay rate. In terms of the notation used in Eq. 22, 
we set 



(fc) lYfik) 
J ij ' 



(23) 



where the index k runs over the sub-components. For 

example, the form factor in B — !■ DtP decays is of the 
form G{z,pjj) = A{z) — p'jjB{z), where pjj, the slope 
of the form factor, is a parameter in the fit and 2; is a 
kinematic variable. The decay rate, which depends on 
the square of G, has terms proportional to 1, pjj and 
{p%)^, thus requiring three sub-components with coef- 
ficients Cj^^ to C^^\ The calculation of the variance 
for the B Dtu component involves the fourth power 
of G and thus requires five sub-components. For the 
B D*£l7 decay we use 18 sub-components to allow the 
fitting of the form factor parameters Ri, R2 and p^, and 
75 sub-components to calculate the associated variance. 
By breaking the components up in this way the fitted pa- 
rameters enter only as multiplicative factors on specific 
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component histograms. M^j , which allows us to use pre- 
made histograms to re-calculate expected yields, avoids 
the need to loop over the simulated events at each step 
in the minimization process and results in a dramatic 
reduction in the required computation time. 



A. Fit parameters and inputs 



TABLE III: Assumed masses, widths, and spins of the four 
hypothetical high-mass states contributing to B — » D^'^TrnW 
decays. 

name mass (GeV/c^) width ( GeV) spin 
ITc 2^61 03 0~ 

X* 2.61 0.3 1 

Yc 2.87 0.1 

Y* 2.87 0.1 1 



TABLE 11: Input parameters for the fit. 



parameter 


value 


Ri 


1.429 ± 0.061 ± 0.044 


R2 


0.827 ± 0.038 ± 0.022 


D" FF slope 


-1.5 ±0.5 


8(3- D^'^T^ev) 


0.0151 ±0.0015 




0.74 ± 0.20 


fo^D-n / DiD^ 


0.87 ± 0.43 




0.36 ±0.30 




0.21 ±0.21 




0.14 ±0.14 


B{B- Z)(*)7r7r^) 


0.011 ±0.011 


B{D*+ D°n+) 


1.7 ±0.4 


0.677 ± 0.005 


B{D° ^ K-n+) 


0.0389 ± 0.0005 


B{D+ K--R+7r+) 


0.0922 ± 0.0021 




1.071 ± 0.009 


/+-//00 


1.065 ± 0.026 



The semileptonic decay widths of i? ^ D£u, B 
D*£U and B — > D**iu are required to be equal for B+ 
and B". We also require isospin invariance in the de- 
cays D** D^*^Tr. As a result, the depend on the 
following quantities: B{B^ D^tU) and form factor 
slope p|, for B D£U and B{B~ D*"£ly) and form- 
factor parameters i?2 and pI, for B ^ D*a7. We 
fix Ri and i?2 to the values obtained in Ref. [9]. The 
background contributions are kept at the values deter- 
mined in the simulation. The overall normalizations of 
the B D'-*K£V and B -> D'-*'^TnTiV components are 
also fixed. For the relevant D decay branching fractions 
we use the values from Ref. [11]. The values of input pa- 
rameters are listed in Table II, where fo* is defined as the 
ratio B{D;+ D°Tr+)/B{D;+ D*°n+) [11, 26], and 

/-I //oo is the ratio of branching fractions B{T{4:S) 

B+B^)/B{T{AS) BW) [11]. All fixed values are var- 
ied in assessing systematic uncertainties. 



B. Fit results 

The fit is performed separately on the electron and 
muon samples. The results of these fits are given in 
Table IV. Both fits give good probabilities. The 
corresponding branching fractions are obtained from 
the B^ results by dividing by the lifetime ratio [11] 



TB-/Tgo = 1.071. The statistical correlations for the 
electron and muon samples are given in Table V. Fig. 3 
and Fig. 4 show the projected distributions on the three 
kinematic variables for the electron and muon samples 
along with the ratio of data over fit. 

The results of the separate fits to the De and _D/i sam- 
ples are combined using the full 8x8 covariance matrix. 
This matrix is built from a block-diagonal statistical co- 
variance matrix, with one 4x4 block coming from the 
fit to each lepton sample, and the full 8x8 systematic 
covariance matrix described in Sec. VI. The systematic 
covariance matrix consists of 4 x 4 matrices for the elec- 
tron and muon parameters and a 4 x 4 set of electron- 
muon covariance terms. The corresponding correlation 
coefficients are given in Table VI. There is an advan- 
tage to combining the electron and muon results after 
the systematic errors have been evaluated; the results are 
weighted optimally {e.g., the difference in lepton identifi- 
cation efficiency uncertainties is taken into account) and 
the from the combination provides a valid measure of 
the compatibility of the electron and muon results. The 
combined results are given in Table IV, and the correla- 
tion coefficients corresponding to the combined statistical 
and systematic errors are given in Table VII. 



C. Fit validation 

The fit was validated in several ways. A large number 
of simulated experiments were generated based on ran- 
dom samples drawn from the histograms used in the fit. 
The fit was performed on these simulated experiments to 
check for biases in the fitted values or associated vari- 
ances. Small biases in the fitted values of several param- 
eters - in no case exceeding 0.1 standard deviations for 
both electron and muon samples - were found. Given the 
smallness of the biases we do not correct the fit results. 
Additional sets of simulated experiments were generated 
with alternative values for the parameters. In each case 
the fit reproduced the alternative values within statistical 
uncertainties. An independent sample of fully- simulated 
events was also used to validate the fit. 

Additional fits were performed on the data to look for 
inconsistencies and quantify the impact of additional con- 
straints. The electron and muon samples were combined 
before fitting; the results were compatible with expec- 
tations. Data samples collected in different years were 
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TABLE IV: Fit results on the electron and muon samples, and their combination. The first error is statistical, the second, 
systematic. 



Parameters 


De sample 


Dfjb sample 


combined result 


pl 

Pl> 

B{D*°eu){%) 


1.23 ±0.05 ±0.08 

1.23 ±0.02 ±0.07 
2.38 ±0.03 ±0.14 
5.45 ±0.03 ±0.22 


1.13 ± 0.07 ± 0.09 1.20 ± 0.04 ± 0.07 
1.24 ± 0.03 ± 0.07 1.22 ± 0.02 ± 0.07 

2.26 ± 0.04 ± 0.16 2.34 ± 0.03 ± 0.13 

5.27 ± 0.04 ± 0.37 5.40 ± 0.02 ± 0.21 


X^/n.d.f. (probability) 


422/470 (0.94) 


494/467 (0.19) 


2.2/4 (0.71) 



TABLE V: Statistical correlation coefficients between param- 
eters from the fits to the electron and muon samples. 





De sample 


Djj. sample 




pI pI, b{d) 


pI pI, b{d) 


pI- 


-0.299 


-0.302 


B{D) 


±0.307 ±0.180 


±0.279 ±0.198 


B{D*) 


-0.388 ±0.075 -0.526 


-0.396 ±0.069 -0.519 



TABLE VI: Correlation coefficients for systematic errors. The 
upper and lower diagonal blocks correspond to electrons and 
muons, respectively. 





De sample 


Dfi sample 




pI pi, B{D) B{D*) 


pI pI, b{d) 




-0.02 




B{D) 


±0.74 ±0.08 




B{D*) 


-0.22 ±0.36 ±0.31 




pI 


±0.75 -0.17 ±0.47 -0.35 




pI' 


-0.07 ±0.98 ±0.02 ±0.31 


-0.15 


B{D) 


±0.46 ±0.00 ±0.64 ±0.17 


±0.16 -0.03 


B{D*) 


-0.17 ±0.19 ±0.12 ±0.54 


-0.48 ±0.17 ±0.67 



TABLE VII: Output correlation matrix for combined samples. 





pI pI, b{d) 


pi- 


-0.129 


rn 


±0.609 ±0.023 


B{D*) 


-0.285 ±0.308 ±0.283 



fitted separately; the fit results agree within statistical 
uncertainties. The minimum number of expected entries 
per bin was varied from 10 to 100; the impact on the 
fitted parameters was negligible. Different binnings in 
the variables p^, pj, and cos 6 b-di were tried; the fit 
results were in each case consistent with the nominal val- 
ues. The boundaries of the D mass peak and sideband 
regions were varied by ±2 MeV/c^; the impact on the 
fitted parameters was negligible. 

Additional fits were performed in which i?i and i?2 
were treated as free parameters. The results, includ- 
ing associated systematic uncertainties, arc given in Ta- 
ble VIII. Correlation coefficients for the combined fit 



are given in Table IX. The three D* form factor pa- 
rameters are highly correlated. Comparing this set of 
parameters with the previous measurement [9], we find 
they are consistent at the 36% C.L. 



VI. SYSTEMATIC STUDIES 

There are several soiirces of systematic uncertainty in 
this analysis. Table X summarizes the systematic uncer- 
tainties on the quantities of interest; these were used in 
determining the systematic errors and correlations given 
in Tables IV and VI. 

The parameters R\ and i?2 are varied taking their cor- 
relation (—0.84) into account. We transform Ri and 
i?2 into a set of parameters R'^ and R'2 that diagonal- 
ize the error matrix, and vary R'l and R'2 independently. 
The D** form factor shape is varied in two ways: the 
slope is varied from —2.0 to —1.0, and the approxima- 
tion Bi from Ref. [23] is replaced with B2 (see also Ap- 
pendix A). The total and relative branching fractions of 
the D** components in _B ^ D^*^Trtr' decays are var- 
ied independently using the values in Table II. The 
D* /D ratio of non-resonant decays, which is defined by 
qNR^B D*wiU)/B'^'^{B DiriU), is 0.3 in the nom- 
inal fit; we vary the ratio from 0.1 to 1.0. The branching 
fraction of S ^ D^*^mriT' decays is varied as given in 
Table II, and the production ratios for the states used 
to model B D(*hTreu decays, X;/Xe, Y^/Y^, X^/Y^ 
and X*/Y*, are varied independently from 0.5 to 2.0. To 
evaluate the effect of Di Dtttt decays [29], one half of 
the B D'-*^T:TTiV component is replaced by Di Dtttt 
decays; the differences in fitted values are taken as sys- 
tematic uncertainties. 

The other parameters listed in Table II are also var- 
ied within their uncertainties. The determination of the 
number of BB events introduces a normalization uncer- 
tainty of 1.1% on the branching fractions. The uncer- 
tainty in the luminosity ratio between on-peak and off- 
peak data is 0.25%. 

The B momentum distribution is determined from the 
well- measured beam energy and B^ mass. The uncer- 
tainty of 0.2 MeV in the beam energy measurement leads 
to a systematic error. Uncertainties arising from the 
simulation of the detector response to charged particle 
reconstruction and particle identification are studied by 
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TABLE VIII: Results on the electron, muon and combined samples when fitting Ri and i?2- 



Parameters De sample Dfi sample combined result 

p5 1.22 ± 0.05 ± 0.10 1.10 ± 0.07 ± 0.10 1.16 ± 0.04 ± 0.08 

Po. 1.34 ± 0.05 ± 0.09 1.33 ± 0.06 ± 0.09 1.33 ± 0.04 ± 0.09 

Ri 1.59 ± 0.09 ± 0.15 1.53 ± 0.10 ± 0.17 1.56 ± 0.07 ± 0.15 

i?2 0.67 ± 0.07 ± 0.10 0.68 ± 0.08 ± 0.10 0.66 ± 0.05 ± 0.09 

B{D°a7){%) 2.38 ± 0.04 ± 0.15 2.25 ± 0.04 ± 0.17 2.32 ± 0.03 ± 0.13 

BlD*°eu){%) 5.50 ± 0.05 ± 0.23 5.34 ± 0.06 ± 0.37 5.48 ± 0.04 ± 0.22 



xVn.d.f. (probability) 416/468 (0.96) 488/464 (0.21) 2.0/6 (0.92) 



TABLE IX: Output correlation coefficients for combined sam- 
ples with Ri and R2 fitted. 





Pd Pd' Ri 


R2 B{D) 


pI, 


-0.435 




Ri 


-0.252 +0.752 




R2 


+0.519 -0.787 -0.740 




B{D) 


+0.602 -0.056 +0.114 


+0.102 


B{D*) 


-0.310 +0.406 +0.139 


-0.309 +0.212 



varying the efficiencies and mis-identification probabili- 
ties based on comparisons between data and simulation 
on dedicated control samples. The uncertainty arising 
from radiative corrections is studied by comparing the re- 
sults using PHOTOS [30] to simulate final state radiation 
(default case) with those obtained with PHOTOS turned 
off. We take 25 % of the difference as an error. The un- 
certainty in the simulation of brcmsstrahlung is based 
on an understanding of the detector material from stud- 
ies of photon conversions and hadronic interactions. The 
uncertainty associated with the charge particle vertex re- 
quirements for the D and B decay points is evaluated by 
loosening the vertex probability cuts. The uncertainty 
arising from B^B'^ mixing is negligible. 

Branching fractions in background simulations are var- 
ied within their measured uncertainties [11]. The inclu- 
sive differential branching fractions versus D momentum 
for B meson decays to D^, D^, £)+ and D~ mesons, 
which affect some background components, are varied us- 
ing the measurements from Ref. [20]. 

The overall covariance matrix for the 8 fitted quan- 
tities (4 electron and 4 muon parameters) is built from 
the individual systematic variations as follows. For each 
variation taken, an 8-component vector Aa of parameter 
differences between the alternative fit and the nominal fit 
is recorded. The ij element of the systematic error co- 
variance matrix is the sum over all systematic variations 
k: 

14,.,, = ^Aaf)Aaf. (24) 

k 

The corresponding correlation matrix is given in Ta- 
ble VI. 



VII. DETERMINATION OF \Vcb\ 

The combined fit results with their full covariance ma- 
trix are used to calculate S(1)|K6| and .F(l)|14(,|: 

g{l)\Vcb\ = (43.1 ± 0.8 ± 2.3) X 10"^ (25) 
J^{l)\Vcb\ = (35.9 ± 0.2 ± 1.2) X 10"^ (26) 

The errors are statistical and systematic, respectively. 
The associated correlations are +0.64 (between ^?(l)|Vc6| 
and p%), +0.56 {H^)\Vcb\ and pi,) and -0.07 {g{l)\Vcb\ 
and T{l)\Vcb\). 

Using the values of .F(l)|146| and 5(1)|T46| given above 
along with calculations of the form factor normalizations 
allows one to determine \Vcb\- Using a recent lattice QCD 
calculation, ^^(1) = 1.074 ± 0.018 ± 0.016 [33], multiplied 
by the electroweak correction [32] of 1.007, we find 

D£iy: \Vcb\ = (39.9±0.8±2.2±0.9) X 10"^ (27) 

where the errors are statistical, systematic and theoret- 
ical, respectively. For B — > D*(17 we use a lattice QCD 
calculation of the form factor, = 0.921 ± 0.013 ± 

0.020 [31], along with the electroweak correction factor, 
to find 

D*eu: \Vcb\ = (38.6 ± 0.2 ± 1.3 ± 1.0) X 10-^(28) 

The fits with Ri and R2 as free parameters give 

g{l)\Vcb\ = (42.8 ± 0.9 ± 2.3) X 10^3 (29) 
T{l)\Vcb\ = (35.6 ± 0.3 ± 1.0) X 10"^ (30) 

with correlation coefficients +0.92 (between C/(l)|T4h| and 
pI), +0.41 {J^{l)\Vcb\ and pi,) and -0.03 {g{l)\Vcb\ and 
m\Vcb\). 



VIII. DISCUSSION 

The branching fractions and slope parameters mea- 
sured here for B DtP and B ^ D* tU are consistent 
with the world averages [25] for these quantities. The 
measurements of pi and ^(l)]!^^] represent significant 
improvements on existing knowledge. The experimental 
technique used here, namely a simultaneous global fit to 
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TABLE X: Systematic uncertainties on fitted parameters, given in %. Numbers are negative when the fitted value decreases 
as input parameter increases. 









Electron sample 










Muon 


sample 






item 


2 

Pd 


2 

Pd* 


B{Deu) B{D*eu) G{1)\Va\ 


^(X)\Vcb\ 


2 

Pd 


2 

Pd' 


B(DW) BiD'iU) g{l)\Vcb\ 




R'l 


0.44 


2.74 


0.71 


-0.38 


0.60 


0.71 


0.50 


2.67 


0.74 


-0.40 


0.63 


0.70 




-0.40 


1.02 


—0.18 


0.30 


—0.32 


0.49 


-0.45 


0.96 


—0.19 


0.30 


—0.33 


0.48 


D** slope 


— 1.42 


—2.52 


—0.07 


—0.09 


—0.82 


—0.87 


— 1.42 


—2.58 


—0.10 


—0.10 


—0.77 


—0.92 


D** FF approximation 


— U.8/ 





—0.12 


0.19 


—0.54 


U.ZU 


nn 
— U.99 


U.o9 


—0.12 


0.21 


—0.59 


Qn 


B{B D'- >n£V) 


0.28 


—0.27 


-0.22 


-0.80 


0.04 


—0.49 


0.59 


— 0.32 


-0.13 


-0.86 


0.24 


—0.54 




—0.39 


0.16 


-0.38 


0.16 


-0.41 


0.13 


—0.50 


0.17 


-0.41 


0.18 


-0.47 


0.15 


fD[D*Ti/DiD* 


-2.30 


1.12 


-1.53 


0.97 


-2.07 


0.85 


-3.13 


1.23 


-1.53 


1.02 


-2.41 


0.93 


1.82 


-1.14 


1.30 


-0.65 


1.65 


-0.70 


2.44 


-1.15 


1.35 


-0.72 


1.91 


-0.75 


Id-k/D* 

fD*Tr/D[ 

NR D*/D ratio 


-0.88 


-1.28 


0.36 


0.17 


-0.31 


-0.34 


-0.83 


-1.23 


0.31 


0.18 


-0.27 


-0.33 


-0.21 


-0.05 


-0.13 


0.21 


-0.18 


0.09 


-0.30 


-0.04 


-0.15 


0.23 


-0.23 


0.10 


0.58 


-0.16 


0.11 


-0.09 


0.38 


-0.04 


0.66 


-0.16 


0.11 


-0.09 


0.40 


-0.03 


B{B- D'-'^TTTveiy) 


1.19 


-1.97 


0.25 


-1.28 


0.78 


-1.28 


1.98 


-1.71 


0.40 


-1.20 


1.20 


-1.18 


X* IX and y lY ratio 


0.61 


-1.15 


0.09 


-0.27 


0.39 


-0.52 


0.74 


-1.02 


0.08 


-0.24 


0.42 


-0.47 


XjY and X*IY* ratio 


0.76 


-0.83 


0.21 


-0.65 


0.52 


-0.60 


1.09 


-0.76 


0.25 


-0.63 


0.68 


-0.57 




2.22 


-1.54 


0.74 


-1.08 


1.63 


-1.05 


2.74 


-1.48 


0.76 


-1.06 


1.81 


-1.03 


B{D*+ PV+) 


-0.14 


-0.01 


-0.10 


0.07 


-0.12 


0.03 


-0.16 


-0.01 


-0.10 


0.07 


-0.13 


0.03 


0.73 


-0.01 


0.43 


-0.34 


0.62 


-0.17 


0.80 


-0.00 


0.41 


-0.33 


0.61 


-0.17 




0.69 


0.02 


-0.21 


-1.63 


0.29 


-0.80 


0.92 


0.12 


-0.27 


-1.68 


0.35 


-0.80 


B{D+ K-7T+-n:+) 


-1.46 


-0.42 


-2.17 


0.30 


-1.89 


0.01 


-1.43 


-0.42 


-2.10 


0.28 


-1.77 


-0.01 


Tb-Itbo 


0.26 


0.16 


0.63 


0.27 


0.46 


0.19 


0.22 


0.16 


0.58 


0.28 


0.41 


0.19 


/+-//00 


0.88 


0.43 


0.66 


-0.53 


0.82 


-0.12 


0.91 


0.48 


0.57 


-0.52 


0.75 


-0.10 


Number of BB events 


0.00 


-0.00 


-1.11 


-1.11 


-0.55 


-0.55 


0.00 


-0.00 


-1.11 


-1.11 


-0.55 


-0.55 


Off-peak Luminosity 


0.05 


0.01 


-0.02 


-0.00 


0.02 


0.00 


0.07 


0.00 


-0.02 


-0.00 


0.02 


-0.00 


B momentum distrib. 


-0.96 


0.63 


1.29 


-0.54 


-1.15 


0.48 


1.30 


-0.10 


1.27 


-0.64 


1.31 


-0.35 


Lepton PID eff 


0.52 


0.16 


1.21 


0.82 


0.90 


0.46 


3.30 


0.06 


5.11 


5.83 


1.99 


2.90 


Lepton mis-ID 


0.03 


0.01 


-0.01 


-0.01 


0.01 


-0.00 


2.65 


0.70 


-0.59 


-0.50 


1.06 


-0.01 


Kaon PID 


0.07 


0.80 


0.28 


0.23 


0.18 


0.38 


1.02 


0.71 


0.35 


0.29 


0.70 


0.39 


Tracking eff 


-1.02 


-0.43 


-3.35 


-2.00 


-2.25 


-1.15 


-0.63 


-0.28 


-3.37 


-2.09 


-2.02 


-1.14 


Radiative corrections 


-3.13 


-1.04 


-2.87 


-0.74 


-3.02 


-0.71 


-0.76 


-0.61 


-0.82 


-0.25 


-0.79 


-0.33 


Bremsstrahlung 


0.07 


0.00 


-0.13 


-0.28 


-0.04 


-0.14 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


Vertexing 


0.83 


-0.64 


0.63 


0.60 


0.78 


0.09 


1.79 


-0.76 


0.97 


0.54 


1.41 


0.01 


Background total 


1.39 


1.12 


0.64 


0.34 


1.07 


0.51 


1.58 


1.09 


0.67 


0.38 


1.16 


0.49 


Total 


6.25 


5.66 


6.01 


4.03 


5.99 


3.20 


8.12 


5.47 


7.35 


7.07 


6.06 


4.23 



B D^Xtu and B D^XtU combinations, is com- 
plementary to previous measurements. In particular, it 
does not rely on the reconstruction of the soft transition 
pion from the £)* — > Dir decay. 



The results obtained here, which arc given in Ta- 
ble IV, can be combined with the existing BaBAR mea- 
surements listed in Table XI. For B D*£p, we com- 
bine the present results with two BaBar measurements 
of and .7-"(l)| Vcf,| [9, 10] and four measurements of 
B{B D*eu)[6, 9, 10]. We neglect the tiny statistical 
correlations among the measurements and treat the sys- 
tematic uncertainties as fully correlated within a given 
category (background, detector modeling, etc.). We as- 
sume the semileptonic decay widths of and B^ to 
be equal and adjust all measurements to the values of 
the T{4:S) and D decay branching fractions used in this 



article to obtain 

B{B- ^ D*°£P) = (5.49 ± 0.19)% (31) 
pI,. = 1.20 ±0.04 (32) 
T{l)\Vcb\ = (34.8 ± 0.8) X 10"^ (33) 

The associated probabilities of the averages are 
0.39, 0.86 and 0.27, respectively. The average of the 
B{B D£U) result with the two existing BaBAR mea- 
surements [6] is 

B{B- -> D"ii7) = (2.32 ± 0.09)% (34) 

with a probability of 0.88. 

The simultaneous measurements of tJ(l)|Vc6| and 
JF(l)|ycfc| allow a determination of the ratio 0(l)/J^(l) 
which can be compared directly with theory. We find 

Measured: g{l)/T{l) =1.20±0.09 (35) 
Theory: g(l)/J^{l) =1.17 ±0.04, (36) 
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TABLE XI: Previously published BABAR results [6, 9, 10]. 



Parameters 


Ref. [9] 


Ref. [10] 


Ref. [6] 


B{B- D*«£I7)(%) 
Z3(B" D*+ft7)(%) 

^(1)1^61 (xlO-3) 


1.191 ±0.048 ±0.028 

4.69 ±0.04 ±0.34 
34.4 ± 0.3 ± 1.1 


1.16 ±0.06 ±0.08 
5.56 ±0.08 ±0.41 

35.9 ±0.6 ±1.4 


5.83 ±0.15 ±0.30 
5.49 ±0.16 ±0.25 


B{B- D"ei>){%) 
B{B° ^D+eu){%) 


2.33 ± 0.09 ± 0.09 
2.21 ±0.11 ±0.12 



where we have assumed the theory errors on J-"(l) [31] 
and ^(1) [33] to be independent. The measured ratio is 
consistent with the predicted ratio. 

The excellent description obtained in this fit, at the 
1% statistical level, of the dominant Cabibbo- favored 
semileptonic decays will facilitate the determination of 
decay rates of Cabibbo-suppressed decays over a larger 
kinematic region than has been feasible to date. This 
will result in a reduction in the theoretical uncertainty 
on the determination of \Vub\- 

To summarize: we use a global fit to D^i and D~^£ 
combinations to measure the form factor parameters 

pIj = 1.20 ±0.04 ±0.07 
pjj. = 1.22 ±0.02 ±0.07, (37) 

in the commonly used HQET-based parameteriza- 
tion [16] and the branching fractions 

B{B-^D°eu) = (2.34 ±0.03 ±0.13)% 
B{B-^D*°iP) = (5.40 ±0.02 ±0.21)%, (38) 

where the first error is statistical and the second sys- 
tematic. The fit assumes the semileptonic decay widths 
of 5+ and B° to be equal. These results are consistent 
with previous BaBAR measurements [6, 9, 10]. From these 
slopes and branching fractions we determine 

g{l)\Vcb\ = (43.1 ± 0.8 ± 2.3) X 10"^ 
J^{l)\Vcb\ = (35.9 ± 0.2 ± 1.2) X 10"^. (39) 

The 5(l)|Vc6| value is twice as precise as the current 
world average. The precision on ^(l)]!^;,] equals that of 
the best single measurement, while coming from a com- 
plementary technique. From these results, we extract two 



values for \Vcb\- 

D*eiy: \Vcb\ = (38.6 ± 0.2 ± 1.3 ± 1.0) X 10"^ 
Diiy: \Vcb\ = (39.9 ± 0.8 ± 2.2 ± 0.9) X 10-^(40) 

where the errors correspond to statistical, systematic and 
theoretical uncertainties, respectively. 
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APPENDIX A: MODELING OF S ^ D**eve 
DECAYS 



The differential decay rates of B D**£p£ decays are 
given as functions of w and 6 [23]. This is the angle 

between the charged Icpton and the charmed meson in 
the rest frame of the virtual W boson. Thus 9 is related 



to Og, which is defined in Fig. 5, such that 
cos^ = cos(7r — 6e) = — cosOg. 



(Al) 



In the following subsections, wo use the same notation 
as above, r and R, for the mass ratios of all four D** 
mesons. However, it is implied that these are the ratios 
taken with corresponding charmed meson masses. The 
following notations are also used in the form factor for- 
mulae in the following subsections : 



1 



£b = 



2mb ' 



£c = 



2mc 



(A2) 



and 



A = energy of the ground state doublet (D and D*) 
A' = energy of the excited |^ doublet {Di and D2 ) 
A* = energy of the excited 5^ doublet {Dq and D[). 

(A3) 



a. B -» DiiU 
The differential decay rate is given by 

= r.r^iw' - 1)1/2 (A4) 

aw a cos 
where = ^'^'^l^'a"" and 

= (1 - cos2 0)[{w - r)fv, + [w^ - 1)(/V3 + rfv,)? 
+(1 - 2rw + r^m + cos^ e){f^^ + {w^ - 1)/1) 

-Acose^/w^-lfvjA\ (A5) 

and fvi{w), fv^i^), fvsiw) and fA{w) are form factors 
which are given by 

V6fA = -{w+l)T 

-eb{w - 1)[{A'+A)t - (2w + 1)ti - T2] 
-SciMwA' - A)t - 3{w - 1)(ti - T2)] 

V6/v, = (1 - w^)t 

-Sb{w^ - 1)[{A' + A)t - {2w + l)ri - t^] 
-ec[4(w + l){wA' - A)t - 3(^2 _ i)(ri - T2)] 

V6fv2 = -3t - 3eb[{A' + A)t - {2w + l)ri - T2] 
-£c[(4w - l)ri -I- 5t2] 

V&fv, = {w- 2)t 

-h£b(2 + w)[{A[_ + A)t - {2w + 1)ti - T2] 
+ec[A{wA' - A)t +{2 + w)ti + {2 + 3«;)r2]. 

,(A6) 

Here t is the leading Isgur- Wise function, which is as- 
sumed to be a linear form [23] 



t{w) = t{1)[1 + t' {w - I)]. 



(A7) 



Uncertainty in first order expansion of Isgur- Wise func- 
tion is parameterized in n and T2. In approximation Bi 

one sets 



Tl=0, T2=0, 



(A8) 
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while in approximation B2 one takes 
Ti = At, T2 = —A't. 

b. s -> D^en 

The differential decay rate is given by 

J2 



(A9) 



^^ = r..3(.._i)3/.ix.,K.), (Alo) 



where 



= ^(1 - cos^ e)[iw - r)kA, + {w^ - l)(fcA3 + rkA,)f 
+(1 - 2rw + r^)[{l + cos^ e){k\^ + {w^ - 1)4) 

-Acose^/w^^ -lkA^kv] (All) 

and kv{w), kAi{w), kA2{w) and kAsiw) are form factors 
which are given by 



kv = -T- SbiiA' + A)t - {2w + 1)ti - T2] 

-Sciri - T2) 
kAi = -{1 + w)t 

-eb{w - 1)[{A' + A)t - {2w + l)ri - T2] 

-ec{w - 1)(ti - T2) 
kA2 = -2ecTi 

kA3=T + eb[{A' + A)t - (2w + l)ri - T2] 
-edn + T2). 



(A12) 



The differential decay rate is given by 



dw d cos 9 



^ = ryiw-" - 1)3/2 J- (A13) 



where 



Jfl. {w, e) = {l- cos2 ^)[(1 + r)g+ - (1 - r)g-f (A14) 
and g+{w) and g-{w) are form factors which are given 

by 



9+ = £c 
9-=C 



wA'—A 



w+1 



2{w - l)Ci - 3C 

A'(2w+l)-Aiw+2)^ _ _ ^^^^ 



w+1 



(A15) 



with 

w + 1 
C(w) = —^t{w). 

In approximation Bi one uses 

Ci = o, C2 = o, 

while in approximation B2 one takes 

Ci = AC,_ C2 = -^*C- 
d. s -> D[eu 

The differential decay rate is given by 
d^Tn' 



(A16) 

(A17) 
(A18) 



dwd cos 6 



^=r,r3(«;2_ 1)1/2 (A19) 



where 

= (1 - cos2 9)[{w - r)gv, + - l)(5y3 + ^ffv^)]' 
+(1 - 2rw + r2)[(l + cos^ e){gl^ + {w' - l)gl) 

-4cos6'^/^t;2 - Igv.gA] (A20) 

and gvi{w), gviiw), gvs,{w) and gA{w) are form factors 
which are given by 



9A = C + e. 

-£b 



uA'-A 



w+1 

A^2w+l)-A(w+2) ^ _ _ 1)^^^ 

gv, ={w-\)Q + e,{wA* - A)C 

-Sb [{A*{2w + 1) - A{w + 2))C - 2(w;2 _ i)^,] 

9V2 = 2£cCi 



5^3 = -C - £ 

' A'{2w+l)-A(w+2) 



wA^—A/- I o/^ 

^+1 C + ^Ci 



+£6 



C-2(m;-1)Ci 



(A21) 



